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We present the first result for the hyperon vector form factor f 1 for S — > T> + W and E~ — > nW 
semileptonic decays from fully dynamical lattice QCD. The calculations are carried out with gauge 
configurations generated by the RBC and UKQCD collaborations with (2+l)-flavors of dynamical 
domain-wall fermions and the Iwasaki gauge action at /3 = 2.13, corresponding to a cutoff a" 1 = 1.73 
GeV. Our results, which are calculated at the lighter three sea quark masses (the lightest pion mass 
down to approximately 330 MeV), show that a sign of the second-order correction of SU(3) breaking 
on the hyperon vector coupling /i(0) is negative. The tendency of the SU(3) breaking correction 
observed in this work disagrees with predictions of both the latest baryon chiral perturbation theory 
result and large N c analysis. 



I. INTRODUCTION 

Greater knowledge of the vector form factor /i in 
AS* = 1 semileptonic hyperon decays paves the way for 
an alternative determination of the element V us of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix in addition 
to kaon semileptonic (-fQa) decays, leptonic decays of 
kaons and pions, and hadronic decays of r leptons [50] , 
A stringent test of CKM unitarity through the first row 
relation \V u d\ 2 + \ V US \ 2 + \V u b\ 2 = 1 can be accomplished 
with the precision of | V us | [2] . A theoretical estimation of 
the vector coupling /i(0) is required to extract V us from 
the experimental rate of hyperon beta decay [33 0] . 

The matrix element for hyperon beta decays, B\ — > 
B 2 W, is composed of the vector and axial-vector transi- 
tions, (B2(p')\V a (x) + A a (x)\Bi(p)) , which are described 
by six form factors: the vector (/1), weak magnetism 
(/2), and induced scalar (fa) form factors for the vector 
current, and the axial- vector (<?i), weak electricity (172), 
and induced pseudo-scalar (#3) form factors for the axial 
current [3]. The experimental decay rate of the hyperon 
beta decay, B x — > B 2 , is given by 



r = 



^(M Bl - M B2 ) 5 (1 3S)\V us \ 2 \f^ B2 (0)\ 2 



x(l + Arc) 



9x 



(0) 



(1) 



where Gp is the Fermi constant measured from the muon 
life time, which already includes some electroweak ra- 
diative corrections [3]. The remaining radiative cor- 
rections to the decay rate are approximately repre- 
sented by A R ,c [5]. Here, M Bl (M B2 ) denotes the rest 
mass of the initial (final) octet baryon state. The el- 
lipsis can be expressed in terms of a power series in 
the small parameter S — {Mb x — Mb 2 ) / {M Bl + Mb 2 ), 
which is regarded as a size of flavor SU(3) breaking [BJ. 
The first linear term in 5, which should be given by 
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-4% 2 (0)5i(0)/ /i(0) 2 ] Bi ^b 2 [31], is safely ignored as 
small as 0(S 2 ) since the nonzero value of the second-class 
form factor g 2 [9] should be induced at first order of the S 
expansion [BJ. The absolute value of 51(0)7/1(0) can be 
determined by measured asymmetries such as electron- 
neutrino correlation [3J [BJ . A theoretical attempt to eval- 
uate SU(3)-breaking corrections on the vector coupling 
/i(0), whose value is given by SU(3) Clebsch-Gordan co- 
efficients in the exact SU(3) limit, is primarily required 
for the precise determination of \V US \. 

The value of /i(0) should be equal to the SU(3) 
Clebsch-Gordan coefficients up to the second order in 
SU(3) breaking, thanks to the Ademollo-Gatto theorem 
(AGT) [TO]. As the mass splittings among octet baryons 
are typically of the order of 10-15%, an expected size of 
the second-order corrections is a few percent level. How- 
ever, either the size or the sign of their corrections is 
somewhat controversial among various theoretical studies 
at present as summarized in Table [Ij A model indepen- 
dent evaluation of SU(3)-breaking corrections is highly 
desired. Although recent quenched lattice studies sug- 
gest that the second-order correction on /i(0) is likely 
negative 0[H], we need further confirmation from (2+1)- 
flavor dynamical lattice QCD near the physical point. 

Our paper is organized as follows. In Sec. II, we first 
summarize the numerical lattice QCD ensembles used for 
this work and then give the details of our Monte Carlo 
simulations. The numerical results are presented in Sec. 
III. We begin with our determination of the scalar form 
factor fs{q 2 ), which will be defined in the later session, 
at finite momentum transfer. We discuss in detail the in- 
terpolation of the form factor to zero momentum transfer 
and also the chiral extrapolation of the hyperon vector 
coupling /i(0). Finally, in Sec. VI, we summarize our 
results and conclusions. 



II. SIMULATION DETAILS 

In this paper, we will present the first result for the hy- 
peron vector form factor for S° — > Yi + lv and £~ — > nW 
semileptonic decays from simulations with 2-1-1 flavors 
of domain wall fermions (DWFs). We use the RBC and 
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TABLE I: Theoretical uncertainties of /i (0) = /i (0) // X SUC3) (0) for various hyperon beta-decays. HBChPT and EOMS-CBChPT 
stand for heavy baryon chiral perturbation theory and covariant baryon chiral perturbation theory with the extended on-mass- 
shell (EOMS) renormalization scheme. 



i ; ; 

Type of result (reference) 


A — > p 


Zj — > n 


a — > A 


S — > Zj 


Bag model [11] 


0.97 


0.97 


0.97 


0.97 


Quark model [12] 


0.987 


0.987 


0.987 


0.987 


Quark model [13] 


0.976 


0.975 


0.976 


0.976 


1/N C expansion [14] 


1.02(2) 


1.04(2) 


1.10(4) 


1.12(5) 


Full 0(p 4 ) HBChPT [15] 


1.027 


1.041 


1.043 


1.009 


Full 0(p 4 ) + partial 0{p 5 ) HBChPT OS] 


1.066(32) 


1.064(6) 


1.053(22) 


1.044(26) 


Full 0(p 4 ) EOMS-CBChPT [T7] 


0.943(21) 


1.028(02) 


0.989(17) 


0.944(16) 


Full 0{p 4 ) EOMS-CBChPT + Decuplet [17] 


1.001(13) 


1.087(42) 


1.040(28) 


1.017(22) 


Quenched lattice QCD [7) [8] 


N/A 


0.988(29) 


N/A 


0.987(19) 



UKQCD collaboration ensembles, which are generated on 
a 24 3 x 64 lattice with two light degenerate quarks and a 
single flavor heavier quark and the Iwasaki gauge action 
at j3 = 2.13 [T5]. The dynamical light and strange quarks 
are described by DWF actions with fifth dimensional ex- 
tent L s = 16 and the domain-wall height of M$ = 1.8, 
which give a residual mass of am ICS » 0.003. Each en- 
semble of configurations uses the same dynamical strange 
quark mass, am s = 0.04, which is close to its physical 
value [18]. We have already published our findings in 
nucleon structure from the same ensembles in three pub- 
lications, Refs. PH2I] . 

The inverse of lattice spacing is a -1 = 1.73(3) 
[a=0. 114(2) fin], which is determined from the fl~ baryon 
mass [18]. Accordingly, the physical spatial extent is ap- 
proximately 2.7 fm, where the nucleon vector form fac- 
tor at low q 2 doesn't suffer much from the finite size 
effect though such effect may influence other nucleon 
form factors [HI [20]. We choose three values for the 
light quark masses, am u ^ — 0.005, 0.01, and 0.02, which 
correspond to about 330 MeV, 420 MeV and 560 MeV 
pion masses [35]. We use 4780, 2350, and 1580 trajec- 
tories separated by 20 trajectories for am u d — 0.005, 
0.01, and 0.02 [TS]. The total number of configurations 
is 240 for am ud = 0.005, 120 for am ud = 0.01, and 80 for 
itTiud — 0.02 as summarized in Table [IT] 

We make four (two) measurements on each configu- 
ration using a single source location, which is located at 
(x, y, z, t) = (6n, 6n, 6n, 16n) with n — 0, 1, 2, 3 (n = 0, 2) 
for am u d = 0.02 (am u d = 0.01 and 0.005), and then they 
are averaged on each configuration in order to reduce pos- 
sible autocorrelations among measurements. The statis- 
tical errors are estimated by the jackknife method on such 
blocked measurements. The quark propagators are calcu- 
lated by gauge-invariant Gaussian smearing at the source 
with smearing parameters (N,w) = (100,7). Details of 
our calculation of the quark propagators are described in 
Ref. 



III. NUMERICAL RESULTS 

A. Scalar form factor fs(q 2 ) at q 2 = i|^ ax 

We focus on vector couplings /i(0) for two different 
hyperon beta-decays, S° — > T, + W and E~ — > nW. These 
decays are simply denoted by S — > E and E — > N here- 
after. We recall that /f^ s (0) = +1 and f^ N (0) = -1 
in the exact SU(3) limit. For convenience in numerical 
calculations, instead of the vector form factor fi{q 2 ), we 
consider the so-called scalar form factor 

. . (2) 

where f$ represents the second-class form factor, which is 
identically zero in the exact SU(3) limit [9]. The renor- 
malized value of fs(q 2 ) at g^ ax = —(M Bl — Mb 2 ) 2 < 
[33] can be precisely evaluated by the double ratio 
method proposed in Ref. [7], where all relevant three- 
point functions are determined at zero three-momentum 
transfer |q| = 0. For the three-point functions, we use 
the sequential source method. We use the source-sink 
separation of 12 lattice units following previous works 
related to nucleon structure (TT)ll21| . Details of the con- 
struction of the three-point functions from the sequential 
quark propagator are described in Ref. |23j . 

Here we note that the absolute value of the renormal- 
ized /s('Zmax) 18 exactly unity in the flavor SU(3) symmet- 
ric limit, where /s(<Zmax) becomes /i(0), for the hyperon 
decays considered here. Thus, the deviation from unity in 
l/s(<Zmax)l is attributed to three types of the SU(3) break- 
ing effect: (1) the recoil correction (g^ ax ^ 0) stemming 
from the mass difference of B\ and £?2 states, (2) the 
presence of the second-class form factor fsiq 2 ), and (3) 
the deviation from unity in the renormalized /i(0). Tak- 
ing the limit of zero four-momentum transfer of fs{q 2 ) 
can separate the third effect from the others, since the 
scalar form factor at q 2 = 0, /g(0), is identical to /i(0). 
Indeed, our main target is to measure the third one. 
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TABLE II: N con {, N sep , and jV mea s denote the number fo gauge configurations, trajectory separation between each measured 
configuration, and the number of measurements on each configuration, respectively. The table contains the pion, kaon, nucleon, 
E-baryon, and H-baryon mass for each ensemble. 



m u d 


iVconf 


iV sep 


■^V mcas 


M n [GeV] 


M K [GeV] 


M N [GeV] 


M E [GeV] 


M a [GeV] 


0.005 


240 


20 


2 


0.3297(7) 


0.5759(8) 


1.140(12) 


1.330(9) 


1.431(6) 


0.01 


120 


20 


2 


0.4200(12) 


0.6064(11) 


1.237(13) 


1.386(12) 


1.465(8) 


0.02 


80 


20 


4 


0.5580(11) 


0.6651(11) 


1.412(10) 


1.501(9) 


1.544(8) 



In Fig. [T] we plot the absolute value of the renormal- 
ized fs(qZnax) as a function of the current insertion time 
slice. Good plateaus are observed in the middle region 
between the source and sink points. The lines represent 
the average value (solid lines) and their 1 standard de- 
viations (dashed lines) over range of 3 < t/a < 8. The 
obtained values of |/s(?max)l) which are naturally renor- 
malized in the double ratio method, as well as g^ax values 
are summarized in Table HTI1 



B. Interpolation to zero four- momentum squared 



The scalar form factor fs{q 2 ) at g 2 > is also calcu- 
lable with nonzero three-momentum transfer (|q| =/= 0). 
To avoid unnecessary repetition, we simply give a refer- 
ence j8], where all the technical details are available. 

We use the four lowest nonzero momenta: q = 2ir/Lx 
(1, 0, 0), (1, 1, 0), (1, 1, 1), and (2, 0, 0), corresponding to a 
q 2 range from about 0.2 to 0.8 GeV 2 . We then can make 
the g 2 interpolation of fs{q 2 ) to q 2 = by the values of 
fsil 2 ) at g 2 > together with the precisely measured 



value of fs(q ) at q 



< from the double ratio. 



In Fig. [2] we plot the absolute value of the renormal- 
ized fs(<r) as a function of q 2 for S — > £ (upper panels) 
and £ -> N (lower panels) at am ud = 0.005 (left), 0.01 
(middle) and 0.02 (right). In this work, we also calculate 
the time-reversal process B2 — > B% as well as B\ — > B2, 
to get more data points in the q 2 > region. Open cir- 
cles are \fs(q 2 )\ at the simulated q 2 . The solid (dashed) 
curve is the fitting result with the seven lowest-g 2 data 
points by using the monopole (quadratic) interpolation 
form [8] , while the open diamond (square) represents the 
interpolated value to q 2 = 0. 

As shown in Fig. [2] two determinations to evaluate 
/s(0) = /i(0) from measured points are indeed consistent 
with each other. Thus, this observation indicates that 
the choice of the interpolation form does not affect the 
interpolated value /i(0) significantly. We simply prefer 
to use the values obtained from the monopole fit in the 
following discussion. 
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FIG. 1: The absolute value of /s° n (<?m ax ) as a function of 
the current insertion time slice. The upper (lower) panel 
is for S — > E (S — > N) decay. In each panel, results for 
am u d = 0.005, 0.01, and 0.02 are plotted from top to bottom. 
The lines represent the average value (solid lines) and their 1 
standard deviations (dashed lines) over range of 3 < t/a < 8. 



C. Chiral extrapolation of /i(0) 

tio of /i(0)//f U( ' 3 ' ) can be parametrized as /i(0) = 
In order to estimate /i(0) at the physical point, we /i(0)//i (0) = 1+ A/, where A/ represents all SU(3)- 
perform the chiral extrapolation of /i(0). The ra- breaking corrections on /i(0). We then introduce the 
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following ratio [71 [5] : 



R A f(M K ,M n ) = 



A/ 



{Ml-Mlf 



(3) 



where the leading symmetry-breaking correction, which 
is predicted by the Ademollo-Gatto theorem, is explic- 
itly factorized out. The remaining dependence related to 
either the higher order corrections of the SU(3) breaking 
or simulated pion and kaon masses is hardly observed 
within the statistical errors as shown in Fig. [3| 

Indeed, if we simply adopt a linear fit form on AR as 
a function of M\ + Ml to extrapolate the value at the 
physical point: 



R Af (M K , M w ) =R + R x • (Af| + M 



(4) 



much as the ratio of Ra / have shown neither any higher- 
order corrections of SU(3) breaking nor the remaining 
M K and M„ dependences. 

In Fig. [4J we present the results of fx (0) (filled circles) 
as a function of the pion mass squared for 2 — > £ (left 
panel) and £ — > N (right panel). In each panel, fitting 
curves indicated by dashed and solid curves represent 
the fitting results with and without the correction for 
the strange quark mass, respectively. The extrapolated 
results of /i(0) at the physical point, which are denoted 
as open circles, agree very well with the AGT fit results 
indicated by filled diamond symbols. Both results are 
tabulated in Table [V] together with the data calculated 
at all simulated quark masses. The statistical errors from 
the AGT fit are rather smaller than those of the global 
fits. 



the resulting coefficient R±, which is approximately zero, 
ensures that the remaining dependence of either Mk or 
M v is negligible at least within the current statistics. 
This observation suggests that although the simulated 
strange quark mass is slightly heavier than the physical 
mass, the corresponding systematic error is likely to be 
small in the chiral extrapolation of R A f- 

We may rather use fitting the data of R A / to a constant 
to estimate the value at the physical point. The two 
fits are mutually consistent, but the latter provides the 
smaller error as shown in Fig. [3] All fitted results are 



also tabulated in Table IV as well as the values of Ra/ 
given at all simulated quark masses. We thus quote the 
value of R A f at the physical point: 



R Af (M^ hys ,M^ s ) 



-0.524(77) for 5 ->• £ 
-0.529(125) for £ -> N, 



(5) 



in (GeV)~ 4 , which is obtained from the latter fit, as our 
best estimate. We evaluate the SU(3)-breaking correc- 
tion A/ via Eq. ^ together with the physical kaon and 
pion masses and then get 



;^ gt (o) = i+a/ 



0.9737(39) for S -> £ 
0.9734(63) forE^V, 



(6) 



which we call the AGT fit result hereafter. Alternatively, 
we may perform a global fit of the data on f\ (0) as mul- 
tiple functions of M\ - Ml and M\ + Ml 



h (0) = C + {Ci + C 2 ■ (Ml + Ml)) ■ {Ml -Ml) 2 , (7) 

whose form is motivated by the AGT fit. Our simula- 
tions are performed with a strange quark mass slightly 
heavier than the physical mass. To take into account 
this slight deviation in this global analysis of the chi- 
ral extrapolation, we simply evaluate a correction using 
the Gell-Mann-Oakes-Renner relation for the pion and 
kaon masses, which corresponds to the quark mass de- 
pendence of pseudo-scalar meson masses at the leading 
order of ChPT. This correction could be accurate in as 



TABLE III: Results for |/F n (<Zmax)l, where gLx 
Mb 2 ) 2 with (Bi,B 2 ) = (H,S) and (E,iV). 


= -{M Bx - 


3 -> E E — > 
mud gLx [GeV 2 ] |/r(9max)l </Lx [GeV 2 ] 


N 

l/r(aLx)l 


0.005 -0.0103(16) 0.9879(71) -0.0360(30) 
0.01 -0.0063(15) 0.9795(55) -0.0223(28) 
0.02 -0.0019(4) 0.9928(16) -0.0080(7) 


1.0166(112) 

1.0108(39) 

1.0013(6) 


TABLE IV: Results for R A f in (GeV) 


-4 


m u d 3 E 


E -» N 


0.005 -0.386(159) 
0.01 -0.706(160) 
0.02 -0.501(114) 


-0.689(281) 
-0.503(229) 
-0.472(193) 


physical point (linear) -0.498(243) 
physical point (average) —0.524(77) 


-0.717(398) 
-0.529(125) 



The excellent agreement observed here between two 
different fitting procedures indicates that the systematic 
uncertainty stemming from the small deviation of the 
strange quark mass appears to be relatively small in the 
AGT fit, where we directly insert the physical kaon and 
pion masses into Eq. ([3| with the weighted average of 
i?A/ in order to determine /i(0) at the physical point. 
However, we conservatively quote the global fit results 
as our final estimates. The differences between two de- 
terminations may be regarded as the reliability of the 
extrapolation to the physical point in our current uncer- 
tainty. Hence our final results are 



A(0) 



+0.9732(66) (7) (5) for S -> S 
-0.9698(106)(15)(36) for S -> N, 



(8) 



where the first error is statistical, and the second and 
third are estimates of the systematic errors due to our 
choice of q 2 interpolation and the reliability of the ex- 
trapolation to the physical point, respectively. Note that 



%.2 0.2 0.4 0.6 0.8 1 %.2 0.2 0.4 0.6 0.8 I % .2 0.2 0.4 0.6 08 i 

q 2 [GeV 2 ] q 2 [GeV 2 ] q 2 [GeV 2 ] 




q 2 [GeV 2 ] q 2 [GeV 2 ] q 2 [GeV 2 ] 



FIG. 2: Interpolation of fs(q 2 ) to q 2 — 0. The upper (lower) panels are for S — > E (E — > TV) decay at am u d = 0.005 (left), 
0.01 (middle), and 0.02 (right). Open circles are |/s(? 2 )| at the simulated q 2 . The solid (dashed) curve is the fitting result by 
using monopole (quadratic) interpolation form, while the open diamond (square) represents the interpolated value to q 2 = 0. 



since we simulate at a single lattice spacing, the system- 
atic error introduced by the lattice discretization is not 
estimated there. 

It is worth emphasizing that the signs of the second- 
order corrections on /i(0) are consistent with what was 
reported in earlier quenched lattice studies [3 H] and pre- 
liminary results from mixed action calculation gl] and 
n f = 2 + 1 dynamical improved Wilson fermion calcula- 
tions [25 . However, we recall that the tendency of the 
SU(3)-breaking correction observed here disagrees with 
predictions of both the latest baryon ChPT result [T7] 
and large N c analysis [TH . 

We additionally remark that the latter has received 
some criticism from Mateu and Pich 4J. They pointed 
out that the large N c fit including second-order SU(3)- 
breaking effects on /i(0) becomes unreliable within the 
present experimental uncertainties |34j . 



IV. SUMMARY 

We have studied the flavor SU(3)-breaking effect on 
hyperon vector coupling /i(0) for the S° —> E + and 
£~ — » n decays in (2+l)-flavor QCD using domain wall 
quarks. We have observed that the second-order correc- 
tion on /i(0) is still negative for both decays at simu- 
lated pion masses of M T = 330 — 558 MeV. The size of 
the second-order corrections observed here is also com- 
parable to what was observed in our DWF calculations 
of Kis decays 27J. Using the best estimate of \V US \ = 



0.2254(6) with imposing CKM unitarity gS], we then 
predict the values \V us fi(0)\ S -+x = 0.2194(8)y us (15)^ 
and |K s /i(0)|s^iV = 0.2186(8) v„, (24) A . The" former 
is barely consistent with a single experimental result of 
|K«/l(0)|s->s = 0.209(27) g5], albeit with its large ex- 
perimental error. However, the latter is slightly devi- 
ated from the currently available experimental result of 
|K s /i(0)|s->;v = 0.2282(49) [3] due to reaching the value 
of /i(0) with an accuracy of less than one percent. We 
plan to extend our research to evaluate the systematic 
uncertainty due to the lattice discretization error and 
also to decrease the reliance on the chiral extrapolation 
using RBC/UKQCD 2+1 flavor DWF dynamical ensem- 
bles at a second, finer, lattice spacing with simulated pion 
masses closer to the physical point. 
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TABLE V: Results for [/i(0)//i w (0)]h->e and [/i (0)//i w (0)]e->jv, where /f uw (0) = +1 for S -> E and /f 1 '(0) = -1 
for E — > N . The final column gives the global fit results for the physical strange mass. 









m ud 




physical point 




q 2 interpolation 


0.005 


0.01 


0.02 


AGT fit 


Global fit 


[/i(0)//f U(3) (0)] s ^ E 


monopole 


0.9808(79) 


0.9742(58) 


0.9914(19) 


0.9737(39) 


0.9732(66) 




quadratic 


0.9806(81) 


0.9744(57) 


0.9910(19) 


0.9730(38) 


0.9731(67) 


[/i(o)//f u(3) (o)] s ^ N 


monopole 


0.9656(140) 


0.9816(84) 


0.9919(33) 


0.9734(63) 


0.9698(106) 




quadratic 


0.9641(140) 


0.9870(74) 


0.9918(31) 


0.9759(57) 


0.9748(99) 




FIG. 3: Chiral extrapolation of -Ra/ for S — > E (left) and E — > N (right). In each panel, solid and dashed lines represent the 
fitting results obtained by a linear and a constant fit, respectively. The dotted vertical lines mark the position of the physical 
point. The open square (circle) symbols, which have been shifted slightly to the left (right), represent the extrapolated values 
at the physical point. In the right panel, quenched DWF results [8] are also included as open diamonds for a comparison. 
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